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Q Abstract. High resolution N-body simulations show that the density pro- 

files of dark matter halos formed in the standard CDM cosmogony can be 
fit accurately by scaling a simple "universal" profile. Regardless of their 
mass, halos are nearly isothermal over a large range in radius, but signifi- 
cantly shallower than r~^ near the center and steeper than r~^ in the outer 
regions. The characteristic overdensity of a halo correlates strongly with 
halo mass in a manner consistent with the mass dependence of the epoch 
of halo formation. Matching the shape of the rotation curves of disk galax- 
1^ ies with this halo structure requires (i) disk mass-to-light ratios to increase 

On systematically with luminosity, (ii) halo circular velocities to be systemat- 

ically lower than the disk rotation speed, and (iii) that the masses of halos 
Qh surrounding bright galaxies depend only weakly on galaxy luminosity. This 

Q offers an attractive explanation for the puzzling lack of correlation between 

|-j luminosity and dynamics in observed samples of binary galaxies and of 

satellite companions of bright spiral galaxies, suggesting that the structure 
of dark matter halos surrounding bright spirals is similar to that of cold 
dark matter halos. 
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1. Introduction 

Dark matter halos are often modelled as isothermal potential wells whose 
depth is usually identified with the velocity dispersion of stars in spheroidal 
galaxies or with the disk rotation speed in spirals. Since velocity disper- 
sion increases with galaxy luminosity, this implies that more massive halos 
should surround brighter galaxies, a hypothesis that, however, does not 
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Figure 1. (a-left) The typical density profiles of CDM halos. The leftmost (rightmost) 
system has a mass M200 = 3 x IQ^^ Mq (M200 = 3 x IO^'^Mq). Arrows indicate the value 
of the gravitational softening in each simulation, (b-right) The characteristic overdensity 
of halos as a function of mass. Masses are expressed in units of the current non linear 
mass corresponding to the standard biased CDM spectrum, M, = 3.3 x 10^'^ Mq. 

seem to be supported by observations. Studies of the dynamics of binary 
galaxies and of spiral/satellite samples have revealed a distinct lack of cor- 
relation between the dynamics and the luminosity of the system (White et 
al. 1983, Zaritsky & White 1994). Motivated by this disagreement, we have 
decided to investigate the structure of dark matter halos in the standard bi- 
ased (5 = 1.6) CDM cosmogony using high-resolution N-body simulations. 
A total of 19 halos with masses ranging from that of a rich galaxy cluster 
(~ IQ^^Mq) to that of a dwarf galaxies (lO^^M©) were selected from large 
cosmological simulations and resimulated individually with higher resolu- 
tion. (I assume Hq = 50 km s~^ Mpc~^ for all physical quantities quoted 
in this paper.) The numerical parameters were carefully chosen so that the 
numerical resolution of each simulated halo was the same, enabling mean- 
ingful comparisons between systems of widely different mass. 

2. Results 

Figure la shows the density profiles of four halos, spanning four orders of 
magnitude in mass. The solid curves are fits using a density profile of the 
form suggested by Navarro, Frenk & White (1995a), 

Pcr^t (r/r,)(l + r/r,)2- ^ > 

Remarkably, all the profiles are very well fit by this simple functional form. 
(Here pcrit is the critical density and is a "scale" radius.) If radii are 
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Figure 2. (a-left) The circular velocity as a function of radius for the 19 halos simulated 
in our series. The curves are truncated at the virial radius, r2oo, of each system, (b-right) 
The result of matching the shape of observed disk rotation curves (dotted lines, from 
Persic & Salucci 1995). The dashed lines are the halo contribution to the total circular 
velocity (solid lines) in each case. Radii are expressed in units of the optical radius of 
each galaxy, defined as 3.2 times the exponential radial scalelength. (M / L)disk values are 
in the 5-band. 

expressed in units of the virial radius, r2007 (which determines the total 
mass of the system, M200) there is a single free parameter in eq.(l); the 
characteristic overdensity Sc- Figure lb shows that this overdensity corre- 
lates strongly with the mass of the halo; less massive systems tend to be 
denser than their larger counterparts. Such a trend is expected in hierar- 
chical clustering scenarios such as CDM, since lower mass scales collapse 
at higher redshift and should therefore have typically higher characteristic 
densities. If we define the formation redshift of a halo of mass M as the first 
time when half of its final mass was in progenitors with individual masses 
exceeding some fraction / of M, we can compute analytically the expected 
dependence of Sc on M by assuming that characteristic overdensity of a 
halo merely refiects the density of the universe at the time of formation, 
i.e. Sc(M) oc (1 -|- ZformiM))'^- The curves in Figure lb show that such 
identification provides a good description of the results of the numerical 
experiments for various values of the parameter /, lending support to the 
conclusion that the characteristic density of a halo is determined primarily 
by its formation redshift. (See Navarro, Frenk & White 1995b for details.) 

The circular velocity, Vc{r) = {GM(r)/ry^'^, of all 19 halos is shown 
in Figure 2a as a function of radius. Consistent with eq.(l), Vc rises near 
the center and declines near the virial radius. Over a wide range in ra- 
dius Vc is almost constant, in agreement with the results of previous, lower 
resolution studies (Frenk et al. 1985, 1988, Quinn et al. 1986). Is the ra- 
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dial behaviour of Vc consistent with the observed rotation curves of galaxy 
disks? To compute the rotation curve of disk galaxies with a given rotation 
speed forming within a CDM halo, only two parameters need to be speci- 
fied; the disk stellar mass-to-light ratio, (M / L)disk, and the halo mass or, 
equivalently, its circular velocity at the virial radius, F200 — GM2qqI ^200- 
(This is because, for a given rotation speed, observations constrain the lu- 
minosity and optical radius of a spiral disk, as well as the typical shape of 
the rotation curve.) Slight modifications to the halo structure caused by 
the presence of the disk can be taken into account by assuming that the 
halo responds adiabatically to the assembly of the galaxy (Barnes & White 
1984, Blumenthal et al. 1986). 

Figure 2b shows the values of these two parameters needed to match 
the typical rotation curves of spiral disks, as given by Persic & Salucci 
(1995). Some trends are clear; faster rotating (brighter) disks require larger 
disk mass-to-light ratios, (M/L)disk i^/^, and the asymptotic halo cir- 
cular velocity is systematically lower than that of the disk. Furthermore, 
bright galaxies, i.e. those with Vrot > 200 km/s, are surrounded by halos 
of very similar mass. Indeed, disks with rotation speeds 200 or 300 km/s 
should be surrounded by halos whose mean circular velocities differ by only 
~ 10% and which can be up to a factor of two lower than the disk rota- 
tion speed. This result agrees well with the estimates of Zaritsky & White 
(1994), who found from their study of the dynamics of satellite galaxies 
that the average circular velocity of the halos of bright spirals is between 
180 and 200 km/s. The weak correlation between halo and disk circular 
velocity shown in Fig. 2b also provides a simple explanation for the lack 
of correlation found between luminosity and dynamics of binary galaxies 
and satellite/primary pairs. This result is especially encouraging, since the 
halo properties we infer were chosen to match the shape of the inner rota- 
tion curves, and do not use any information about dynamics at the much 
larger radii probed by binaries or satellite companions. We conclude that 
the structure of dark halos surrounding spiral galaxies is quite similar to 
that of Cold Dark Matter halos. 
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